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gURSTI J UTE SPECIFICATION 

PATENT 

INSTITUT FRANgAIS DU PETROLE 
METHOD AND SYSTEM INTENDED FOR REAL-TIME ESTIMATION OF THE 
FLOW MODE OF A MULTIPHASE FLUID STREAM 
5 AT ALL POINTS OF A PIPE 

inventors : Emmanuel DURET, Eric HEINTZE and Isabelle REY-FABRET 

ABSTRACT 

10 The Invention is a method and system for real-time estimation of the flow 

mode, at all points of a pipe whose structure ss defined by a certain number of 
structure parameters, of a multiphase fluid stream defined by several physical 
quantities providing simplified implementation of hydnDdynamic modules that can be 
integrated In modelling tools. A non-linear neural network is formed with an input 

15 layer having as many inputs as there are stmcture parameters and physical 
quantities, an output layer with as many outputs as there are quantities necessary for 
estimation of the flow mode and at least one Intennediate layer. A learning base Is 
created with predetermined tables connecting various values obtained for the output 
data to the corresponding values of the input data, with iterative determination of the 

20 weighting factors of the activation function allowing to properly connect the values in 
the input and output data tables. In order to avoid singularities of the network output 
data likely to distort the determination of the weighting factors, a sorting procedure is 
used to eliminate non-pertinent data. The main advantages of the method are: 
modelling simplification and time saving. 
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BACKGROUND OF THE INVENTION 

Field of the invention 

The present invention relates to a method and to a system intended for real- 
time estimation of the flow mode of a multiphase fluid stream at all points of a pipe, 
comprising using neural networks. 



Description of the Prior Art 

Transporting hydrocarbons fn3m production sites to treating plants constitutes 
an Important link in the petroleum chain. It is a delicate link because of the complex 

10 interactions between the phases forming the transported effluents. The basic 
objective for operators is to reach an optimum productivity under the best safety 
conditions. The operator therefore have to control as best they can the velocity and 
the temperature so as to avoid unnecessary pressure drops, unwanted deposits and 
unsteady flows. The method that is generally used is to model in the best possible 

15 way the transportation of complex multiphase streams so as to provide at all times 
an image of the flows in the various parts of the production chain, by taking into 
account the precise constitution of the efHuent. the flow rates, the pressures and the 
flow modes. 

There are currently various software tools for simulating the transport of 
20 complex multiphase streams, allowing the design of suitable production equipments 
at an early stage. 

U. S. Patent 5,550.761, French Patent 2,756,044 (U. S. Patent 6,028,992) 
and French Patent 2,756,045 (U. S. Patent 5»960,187) filed by the assignee notably 
describe modelling methods and tools allowing simulation of the transport of complex 
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multiphase streams for steady or transient flow and capable of accounting for 
instability phenomena that occur because of the inregular geometry of the fonnation 
crossed by the pipe or of the topography thereof, which is refen^ed to by specialists 
as "terrain slugging" or "severe slugging". 
5 The simulation tools are as complex as the modelled phenomena. Precision 

and performance can only be obtained after a relatively long modelling time, which is 
not really compatible with real-time management. That is the reason why these 
modelling tools cannot be used as they are for real-time management of the 
production. It therefore appears necessary to use modelling methods offering a good 
10 compromise between calculating speed and.accuracy of results. 

Neural networlcs define a data processing mode simulating the functioning of 
biological neuron systems. In such networks, an element carries out a relatively 
simple calculation such as a weighted sum of the signals present at Its inputs applied 
to a non-linear function, which determines the state of its output- A large number of 
15 such elements interconnected in series and in parallel is used. Proper selection of 
the weightir>g factors allows the network to canry out complex functtons. Networks 
known as retropropagatron networks for example use multiple layers of elements 
defined above. Adaptation of such a network to a precise task is done by "training" 
the network with a certain number of examples and by adjusting the weighting 
20 factors for each element to the suitable values. Input values are presented to the 
network, the output value produced by the network is analyzed and the weighting 
factors are modified to best minimize the difference between the effective value at 
the output and the expected value in the selected example. After a sufficient training 
period, the network is suited to respond to new input values for which the output 
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value is not known a priori and to produce a suitable output value. In its principle, a 
neurai network wori<3 according to a non-linear regression method which is all the 
more effective in relation to conventional methods. 

Such networks are used in many fields such as image recognition, solution of 

5 optimization problems, etc. Initiaiiy, the neural networic is a method suited for 
automatic classification, hence its use In particular for pattern recognition. For these 
applications, two types of networks can be used, the MLP (Muiti Layer Perceptron) 
or the Kohonen networks, well-known to specialists. 

The prior art in the field of neural networics is Illustrated by the following 

10 references: 

Dreyfus G., «Les reseaux de neurones»; M^canlque Industrielle et Mat6riaux, 
n 51, Sept. 98, 

Lippman R.P., An Introduction to Computing with Neural Nets; IEEE ASSP 
IVlagazine, April 1987, or 
15 Pinkus A., Approximation Theory of the MLP Model in Neural Networks; Acta 

Numerica 1999. 

Networics are now also used for non-linear modelling of static data or of 
dynamic processes. The MLP networks are mostly used in this case. This approach 
cun-entiy concerns fields of application such as, for example, anomaly detection or 
20 stock-exchange prediction. 

An example of use of neural networks is described for example in French 
Patent 2,786,568 filed by the assignee. 
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SUMMARY OF THE INVENTION 



In the context of production management, the phenomena to be modelled in 
real time are highly non-linear and the parameters Involved in the modelling 
procedure are numerous. 
5 The invention is a method allowing, alone or in parallel with the 

aforementioned modelling methods, real-time management of the parameters of a 
fluid circulation by using neural networks. 

The invention Is a method for real-time estimation of the flow mode, at all 
points of a pipe whose structure can be defined by a certain number of structure 
10 parameters, of a multiphase fluid stream defined by several physical quantities and 
comprising at least a liquid phase and at least a gas phase, whereas modelling of 
the flow mode comprises: 

fonning a non-linear neural network with an Input layer having as many inputs 
as there are structure parameters and physical quantities, an output layer iiaving as 
15 many outputs as there are quantities necessary for estimation of the flow mode and 
at least one intermediate layer; 

creating a learning base with predetennined tables connecting various values 
obtained for the output data to the con^spondlng values of the input data; and 



20 proper connection of the values in the input and output data tables. 

The method also preferably comprises an analysis of the output data of the 
neural network allowing to sort, among the values of the output data of the neural 
networi^, only the pertinent data to be taken Into account for iterative determination of 
the weighting coefficients of the activation function. 



detennining by iterations weighting factors of the activation function allowing 
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A neural network expresses a continuous mathematical funcUon between 
Input and output data. As it continuously produces results through direct calculations 
that permanently connect Input and output data, the discontinuities linked with the 
solution of systems of equations are consequently avoided, hence an appreciable 
5 time saving and simplified results. The non-linear nature of the activation functions of 
neural networks makes them perfectly suited for modelling of non-linear systems. 

The method comprises for example fomiing a totally connected network and 
using linear output neurons. 

An identity function Is for example selected as the activation function. 
10 The system according to the invention allows real-time estimation of the flow 

mode, at all points of a pipe whose structure can be defined by a certain number of 
structure parameters, of a multiphase fluid stream defined by several physical 
quantities and comprising at least a liquid phase and at least a gas phase. It 
comprises: 

15 means for determining characteristics of a non-linear neural network with an 

input layer having as many inputs as there are stmcture parameters and physical 
quantities, an output layer having as many outputs as there are quantities necessary 
for estimation of the flow mode and at least one intermediate layer. 

means for storing a learning base with predetermined tables connecting 

20 various values obtained for the output data to the corresponding values of the input 
data; and 

means for determining by iterations weighting factors of an activation function 
allowing proper connection of the values in the input and output data tables. 
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The system preferably comprises means for analyzing output data of the 
neural network allowing to sort, among the values of the output data of the neural 
network, only the pertinent data to be taken into account for Iterative detenmination of 
the weighting coefficients of the activation function, 
5 The system of this technique can be used for example in place of a tool 

developed with a conventional technique (and therefore relatively slow) in order, to 
provide results compatible with real-time stream management. The temns of this 
implementation are developed hereafter, 

-10 BRIEF DESCRIPTION OF THE DRAWINGS 

Other features of the method according to the invention will be dear from 
reading the description hereafter, with reference to the accompanying drawings 
wherein: 

Fig. 1 shows an example of a neural network formed to connect 
15 hydrodynamic input data to output data relative to flow modes, and 

Figs. 2 and 3 illustrate two different topological neighbourhood cases leading 
to a different qualification of output data in a validation process of the neural network. 

DETAILED DESCRIPTION OF THE INVENTION 

20 1) Context 

A circulation of multiphase fluids is considered in a pipe vAth at least a liquid 
phase and at least a gas phase, and forming a neural network is attempted allowing, 
from a certain number of geometrical and physical input data relative to the pipe and 
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of physical data relative to the fluids, to give instantly, for each section of the fluid 
stream, an estimation of the flow mode, 

2) inpvt and output data 

The input data are for example: 
5 geometrical data of the network: diameter roughness and angle of inclination 

of the pipe, 

input data qualifying the effluents: density of the gas, density of the liquid, 
viscosity of the gas, viscosity of the liquid, etc.; 

Input data characterizing the mixture: gas/liquid surface tension, volume 
1 0 fraction of gas, barycentric velocity of the mixture. 

The data that the network will evaluate and deliver at two main outputs are: 
dV, the velocity difference between gas and liquid, and 
p, the stratified flow fraction in the pipe section where the flow type is to be 
detemnined; P g [0 ;1]. 

15 Other quantities qualrfying the flow type can be calculated from these two 

outputs. 

3) Structure of the network 

In order to connect the input data to the output data, a preferably IWLP type 
neural network, well-known to the man skilled in the art, is formed since it is 
20 particularly well-suited for physical phenomena modelling. In fact, its structure allows 
describing the dynamic as well as the static components of the phenomena, even by 
fixing, if necessary, some of its parameters at a reified value, therefore physically 
representative. Thus, in the example, knowing physical equations that govern the 
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flows allows enriching the network and to best adapt it to the physical phenomena 
modelled thereby. 

The neural network comprises (Fig,1) three layers : the input layer having ten 
neurons corresponding to the ten data (mentioned above) of the complete physical 
5 model, an output layer having two neurons conresponding to the two parameters dV 
and p sought, and an intemriediate layer, refferred to as hidden layer, whose number 
of neurons No Is optimized. The network is totally connected. The non-linearity of this 
network is obtained by a sigmoid activation function governing the behavior of the 
neurons in the hidden layer. The neurons of the output layer can be selected linear 
10 or non-linear. The activation function can be the Identity function for example. 
4^ Learning: principle and implementation in the example 

a) Principle 

The weights of this structure are determined at the end of a leaming stage; 
this stage supplies the network with a set of data forming the leaming base of the 

15 network, and in optimizing the weights of the network by minimizing the errors noted 
for all the samples of the base, between the output data resulting from network 
calculations and the data expected at the output, given by the base. The errors can 
be the absolute errors between the input and output quantities or the relative errors, 
according to the performance desired for the network. 

20 The generalization powers of the network are then tested from Its capacity to 

properly calculate the two outputs for inputs that are unknown thereto. 

b) Implementation 

!n practice, t>esides the difficulty in selecting the various elements making up 
the networi^, the implementation of the network requires an extensive analysis of the 
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data that make up the teaming base. In fact, even if a network is properly elaborated 
for a given problem, it can give bad results insofar as the learning base that it is 
supplied with contains elements that disturb the optimization of all its weights. !n the 
particular case of fluid flows, the problem is cmcial; in fact, flow mode calculations 

5 carried out from complete and precise models generate a highly Inhomogeneous 
database: for example, a result indicating a stratified flow can be "drowned" in 
neighboring points representing a dispersed flow: this point therefore corresponds to 
either a physical phenomenon that the neural network must be able to find or a 
"singleton" that the network has to smooth on the contrary in order to ensure the 

10 continuity of the result. The difficulty of prior data analysis thus lies in the 
determination of the nature of such points and, on a larger scale, of the various 
points of the base. 

A method intended for a posteriori analysis of the network has therefore been 
elaborated, which identifies and distinguishes the available data representing the 

15 particular physical phenomena to be modelled from the calculation results to be 
ignored. This method produces "compacts" around each singular point considered in 
the n-dimensional space of the inputs (n = 10 in the example described) in order to 
evaluate the degree of isolation of their behavior, thus, the singular points and the 
non-singular points are counted for each "ball" created, the proportion of singular 

20 points in the compact giving then a measure of their isolation. The terms "compact" 
and "ball" are taken in the topological sense of the word. 

Thus, point SI in Fig. 2, which is at the center of a ball and separated from 
the nearest singular point S2 by a neighborhood of many regular points Mi, can be 
considered to be a singleton and thus regularized. On the other hand, the immediate 
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neighborhood of point SI in Fig. 3 has one or more singular points S2. S3» Sk. In 
such a case, one must consider that points SI to Sk are the expression of a 
particular physical phenomenon that has to be taken into account by means of an 
adaptation of the neural network. 
5 This selection method allows identification of the necessary points for 

Infomiing the network sufficiently during the learning stage. The suitable point base 
thus formed can then give infonmation thereto so as to best optimize its parameters. 
5) Results 

The tools and the methods described above allow obtaining; 
10 a database cleared of any point with an abnormal behavior and containing 

only the infonmation-rich points so as to best define the physical system to be 
identified, 

a network capable of interpolating the various flow types encountered so as to 
no longer generate any discontinuity, and also capable of taking Into account the 
15 particular physical phenomena to be modelled despite their atypical behavior, and 
finally 

a tool giving a real-time estimation of the main hydrodynamic information. 
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PATENT 

INSTITUT FRANCAIS DU PETROLE 
METHOD AND SYSTEM INTENDED FOR REAL-TIME ESTIMATION OF 
THE FLOW MODE OF A MULTIPHASE FLUID STREAM 
5 ATALLFOBSTSOFAPIPE 

Inventors : Emmanuel DURET, Eric HEINTZE and IsabeUe REY-FABRET 

ABSTRACT 

to Metho dThe invention ia a method and system intoudod f or real-time estimation of 

the flow mode, at all points of a pipe whose stmcttue is defined by a certain number of 
stmcture parameters, of a multiphase fluid stream defined by several plrysical quantities 
y^Tnvi'Hi'Ti fy fiimpltfied imr?lementation of h vdrodvnamic modules that can bO int<^gr»ted in 
modelling tools . 

1 3 A non-linear neural network is formed with an input layer having as many inputs 

as there are structure parameters and physical quantities, an output layer with as many* 
outputs as there arc quantities necessary for estimation of the flow mode and at least one 
intermediate layer. A learning base is created with predetermined tables connecting 
various values' obtained for the output data to the corresponding values of the input data, 

20 with iterative determination of the weighting factors of the activation ftmction allowing 
to properly connect the values in the input and output data tables. In order to avoid 
singularities of the network output data likely to distort the determination of the 
weighting factors, a sorting procedure is used to eliminate non-^)ertinent data. The main 
advantages of the method are: modelling simplification and time saving. 

25 - Applications : simplified implcmontation of hydrodynam a o — modulOG that can - b e 
intogmttad in mod e llin^g - teolo for example. 
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RAnKGROlJND OF THE INVEim QN | 

Field rtf the Invention 

The present invention lelates to a method and to a system intended for real-time 
estimation of the flow mode of a multiphase fluid stream at all points of a pipe. 
5 comprising using neural networks. 

DACKGROm^JD OF THE DJ\rET^JTlOM 
Description of the Prior Art 

Transporting hydrocarbons from production sites to treating plants constitutes an 
important link in the petroleum chain- It is a delicate link because of the complex 

10 interactions between the phases forming the transported effluents. The basic objective for 
operators is to reach an optimum productivity under the best safety conditions. Tfeeylhe 
opeiutor therefore have to control as best they can the velocity and the temperature so as 
to avoid mmecessary pressure diops» unwanted deposits and unsteady flows. The method 
that is generally used conr i intr. in modeilaR gis to model in the best possible way the 

15 transportation of complex multiphase streams so as to provide at all times an image of 
the flows in the various parts of the production chain, by taking into account the precise 
constitution of the effluent, the flow rates, the pressures and the flow modes. 

There are currently various software tools for simulating the transport of complex 
multiphase streams^ allowing tedie design_of suitable production equipments at an early 

20 stage. 

Patonto US U . S. Patent 5,550,76 K FRr -French Patent_2 ,756.044 (US -U. S. Patent 
6,028,992) and ^French Patent 2,756,045 f^S -U. S. Patent 5,960.187) filed by the 
applican t assigDee notably describe modelling methods and tools allowing ^ 
Si mula Ic s imulation of the transport of complex multiphase streams e<»for steady or 
25 transient flow and capable of taking mto aoooun ta ccounting for instability phenomena 
that occur because of the irregular geometry of the formation crossed by the pipe or of 
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the topography thereof, which is re ferred to by specialists as <tenain slugging^l or 
M'^^€ivere sliigging!!.. 

The simulation tools are as complex as the modelled phenomena. Precision and 
pprfnrmunQQn performance can only be obtained after a relatively long modelling time. 

5 which is not really compatible with real-time management That is the reason why these 
modeUing tools cannot be used as they are for real-time management of Ae production. 
It therefore appears necessary to use modelling methods offering a good compromise 
between calculating speed and accuracy of results, 

¥lio main objoctivc of tho invention io a method allowing, alono or in parall e l 

10 with the qforcmontioncd modelling mcthodg, real tim e managomcnt of tho paramet e rs of 
a fluid oirculation by ujing neural nanvork ft? 

Tt oan b e reminded that neuralN euial networks define a data processing mode 
simulating the functioning of biological neuron systems. In such networks, an element 
carries out a relatively simple calculation such as a weighted sum of the signals present 

15 at its inputs applied to a non-linear function, which determines the state of its output A 
large number of such elements irrterconnected in series and in parallel is used. Proper 
selection of the weighting factors allows the network to carry out complex functions. 
Networks known as retropropagation networks for example use multiple layers of 
^elements defined above. Adaptation of such a network to a precise task is done by 

20 ^"trainin^" the network with a certain number of examples and by adjusting the 
weighting factors for each element to the suitable values. Input values arc presented to 
the network, the output value produced by the network is analvse d analvzed and the 
weighting factors are modified to best minimize the difference between the effective 
value at the output and the expected value in the selected example. After a sufficierrt 

25 training period* the network is suited to respond to new iiq>ut values for which the output 
value is not known a priori and to produce a suitable output value. In its principle, a 
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neural netwoik works according to a non-linear regression method which is all the moie 
effective in relation to conventional methods. 

Such netwodcs are used in many fields such as image recognition, solution of 
opthnization problems, etc. Initially, the neural netwoik is a method suited for automatic 
5 classification, hence its use in particular for pattern recogmtion. For these appHcations, 
two types of networks can be used, the MLP (Multi Layer Perceptron) or the Kohonen 
networks, well-known to specialists. 

The prior art in the field of neural networks is illustrated by the following 

lefbrences: 

10 -Dreyfus G.» <<Lcs r^seaux de neurones»; Mecanique Industrielle et Materiaux, n' 

51,septSeet_98, 

-Lippman R.P., An Introduction to Computing with Neural Nets; IEEE ASSP 
Magazine, April 1987, or 

-Pinkus A., Approximation Theory of the MLP Model in Neural Networks; Acta 
15 Numerica 1999. 

Networks are now also used for non-linear modelling of static data or of dynamic 
processes. The MLP networks are mostly used in this case. This approach currently 
concerns fields of application such as, for example, anomaly detection or stock-exchange 
prediction. 

20 An example of use of neuzal networks is described for example in potent FR A 

French Patent 2 .786.568 filed by the applican t assignee . 

SUMMARY OF THE INVENTION 
In the context of production management, the phenomena to be modelled in real 
time are highly non-linear and the parameters involved in the modelling procedure are 
25 numerous. 
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Th^ ir.v^tinn is a ynftthod allowbi p, alone or in naTallel with the afoicincnt^pa^4 
modelling meibr^^^ real-time maimgement of the pataTTiet^ of a fluid circul^tiPA by 

u^infi neural networl^^ 

The ol:ijoct of the invention is a method »iteadea^for real-time estimation of the 
5 flow mode, at all points of a pipe whose structure can be defined by a certain number of 

structure parameters, of a multiphase fluid stream defined by several physical quantities 

and comprising at least a Kquid phase and at least a gas phase, choractDrized m 

tha twhereas modelling of the flow mode comprises: 

-forming a non-linear neural network with an input layer having as many inputs 
10 as there are structure parameters and physical quantities, an output layer havmg as many 

outputs as there are quantities necessary for estimation of the flow mode and at least one 

intermediate layer?! 

-creating a learning base with predetermined tables connecting various values 
obtained for the output data to the corresponding values of the mpnt data^; and 
15 -determining by iterations wei^ting factors of the activation function allowing *o 

proporly oomi e o tp rODer connection of the values in the h^ut and output data tables. 

The method also preferably comprises an analysis of the output data of the neural 
networic allowing to sort, among the values of the output data of the neural network, only 
the x>ertinent data to be taken into account for iterative determination of the weiglhting 
20 coefficients of the activation function. 

A neural network expresses a continuous mathematical function between input 
and output data. As it continuously produces results through direct calculations that 
permanently coimcot input and output data, the discontinuities linked with the solution of 
systems of equations are consequently avoided, hence an appreciable time saving and 
25 simplified results. The non-linear nature of the activation functions of neural networks 
makes them perfectiy suited for modelling of non-linear systems. 
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The method comprises for example fonning a totaUy comiectsd network and 
using linear output neurons. 

An identity function is for example selected as the activation function. 
The system according to the invention allows leal-time estimation of the flow 
5 mode, at all points of a pipe whose stractuie can be defined by a certain number of 
structure parameters, of a multiphase fluid stream defined by several physical quantities 
and (5ompriBing at least a liquid phase and at least a gas phase. It comprises: 

-means for determining *he-characteristics of a non-linear neural network with an 
input layer having as many inputs as there are structure parameters and physical 
10 quantities, an output layer having as many ou^uts as there are quantities necessary for 
estimation of the flow mode and at least one intermediate layer?; 

-means for storing a learning base with predetermined tables connecting various 
values obtained for the output data to the corresponding values of the input data^ and 

-means for determining by iterations weighting factors of an activation function 
15 allowing to properly - fe^mnoc t ptoper connection of the values in the input and output data 
tables. 

The system preferably comprises means for nnni-y-nin^ nnai Y yin ^ output data of the 
neural network allowing to sort, among the values of the output data of the neural 
network, only the pertinent data to be taken into account for iterative determination of 
20 the weighting coefScients of the activation function. 

The system ^voloppod by of this technique can be used for example in place of a 
tool developped with a conventional technique (and therefore relatively slow) in order to 
provide results compatible with real-time stream management The terms of this 
implementation are dovoloppod develoned hereafter. 

25 
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BRIEF DESCRIPTTON OF THE DRAWINGS 
Other features of the method according to the invention will be clear fion» 
leading the description hereafter, with reference to the accompanying drawings wherein: 
Figur eEis^ 1 shows an example of a neural network formed to connect 
5 hydrodynamic input data to output data relative to flow modes, and 

Figure s Figs. 2 and 3 illustrate two different topological neighbourhood cases 
leading to a drfferent qualification of output data in a validation process of the neural 
network. 

10 DETAILED DESCRIPTION QF THE INVENllQN j 

1) Context 

We conoid e r aA circulation of multiphase fluids is considered in a pipe with at 
least a liquid phase and at least a gas phase^ and we try to tbrm forming a neural network 
is attempted allowing, &om a certain nmnber of geometiica] and physical input data 
15 relative to the pipe and of physical data relative to the fluids^ to give instantiy, for each 
section of the fluid stream, an estimation of the flow mode. 

2) hiput and output data 

The input data are for example: 

-geometrical data of the network: diameter^ roughness and angle of inclination of 
20 the pipe, 

-input data qualifying the effluents: density of the gas, density of the liquid, 
viscosity of the gas, viscosity of the Hquid, etc.; 

-input data characterizing the mixture: gas/liquid surface tension, volume fraction 
of gas, barycentric velocity of the mixture, 
25 The data that the network will evahiate and deliver at two main outputs are: 

-6V^ the velocity difference between gas and liquid, and [ 
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-p, the Stratified flow fraction in flie pipe section where the flow type is to be 
detennincd; 3 e [0 

Other quantities quaKfying the How type can be calculated from these two 
outputs. 

5 3) Structure of the network 

In order to connect the input data to the output data, a preferably MLP type 
neural networlc, weU-known to tfac man skilled in the ait. is formed since it is particularly 
well-suited for physical phenomena modellmg. In feet, its structure allows te 
desGHbedescriimg the dynamic as well as the static components of the phenomena, even 

10 by fixing, if necessary, some of its parameters at a reified value, therefore physically 
representative. Thus, in the example, knowing physical equations diat govern tiie flows 
allows te-eaH6h enriching the network and to best adapt it to the physical phenomena 
modelled thcreby. 

The neural network conq>rise5 (Fig.l) three layers : the input layer oonsisiiD g 
15 e fliaving ten neurons corresponding to fte ten data (mentioned above) of the complete 
physical model, an output layer e o as i s tmp o f having two neurons corresponding to the 
two parameters dV and p sought, and an intermediate layer, referred to as hidden layer, 
whose number of neurons Nc is optimized. The network is totally cormected. The non- 
linearity of this network is obtained by a sigmoid activation fimction governing the 
20 b ohaviom -b ehavior of the neurons in the hidden layer. The neurons of the output layer 
can be selected linear or nou-linear. The activation function can be the identity function 
for example. 

4) Learning: principle and implementation in the example 
a) Principle 

25 The weigihts of this structure are determined at the end of a learning stage; this 

stage coDfiiGts fn stfppjy mgsupplies the network with a set of data forming the learning 
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base of the network, and in optiimzmg the weights of the netwotk by miiiimizmg the 
errors noted for aU the samples of the base, between ttie output data resulting from 
network calculations and the data e;q)ected at the output, given by the base. The errors 
can be the absolute errors between the input and output quantities or the relative errors, 
5 according to the performance desired for the netwoik. 

The generalization powers of the network are then tested from its capacity to 
properly calculate the two outputs for inputs that are unknown thereto, 
b) Implementation 

In practice, besides the difficulty in selecting the various elements making up the 

10 network, the implementation of the network requires an extensive analysis of the data 
that make up the learning base. In feet, even if a network is properly elaborated for a 
given problem, it can give in fine b ad results insofar as the learning base that it is j 
supplied with contains elements fliat disturb the optimization of all its weights. In the 
particular case of fluid flows, the problem is crucial; in feet, flow mode calculations 

15 carried out from complete and precise models generate a highly inhomogcneous 
database: for example, a result indicating a stratified flow can be «^drowned»;i in 
noighbouri ngneighboring points representing a dispersed flow: this point therefore 
corresponds to either a physical phenomenon that the neural network must be able to find 
or a -e^lsingleton^ll. that the network has to smooth on the contrary in order to ensure the 

20 continuity of the result The difSculty of prior data analysis thus lies in the determination 
of the nature of such points and, on a larger scale, of the various points of the base. 

A method intended for a posteriori analysis of the network has therefore been 
elaborated, which identifies and distinguishes the available data representing the 
particular physical phenomena to be modelled from the calculation results to be ignored. 

25 This method ooriiaists in producin g pfoduces <'compacts»" around each singular point 
considered in the n-dimensional space of the inputs (n = 10 in the example described) in 
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order to evaluate the degree of isolation of their behaviour, thxis, the singular points and 
the non-singular points are ooxinted for each 'W^all^^" created, the proportion of singular 
points in the compact giving then a measure of their isolation. The terms ^<lconipac^;;i 
and «!lball»/ are taken in the topological sense of the word. 

5 Thus, point SI in Fig._2, which is at the eeateecenter of a ball and separated from 

the nearest singular point S2 by a neighbourho ^neighborhood of many regular points 
Mi, can be considered to be a singleton and thiis regularized. On the other hand, the 
immediate noighbourh oed neighborhood of point SI in Fig._3 e onaioto of has one or more 
singular points S2, S3, Sk. In such a case, one must consider that pomts SI to Sk are 

10 the e«pression of a particular physical phenomenon that has to be taken into account by 
means of an adaptation of the neiural network. 

This selection method allows lo idcntif y identification of the necessary points for 
informing the network sufficiently during the learning stage. The suitable point base thus 
formed can then give information thereto so as to best optimi2:e its parameters. 

15 5) Results 

The tools and the methods described above allow to obtai n obtaining : 
-a database cleared of any point with an abnormal bohoviou rb ehavior and 
containing only the iufotmation-rich points so as to best define the physical system to be 
identified, 

20 -a network capable of interpolating the various flow types encomitered so as to nQ 

longer g enerate aeanv discontinuit y an y l ong e r , and also enable of taking into accoimt 
the particular physical phenomena to be modelled despite their atypical 
behQvioui b Qhavior. and finally | 
a tool giving a reaVthne estimation of the main hydrodynamic infomiation. 
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